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^ , Abstract 

' ' Relativistic jets are ubiquitous in astrophysical systems that con- 

;_ I ■ tain compact objects. They transport large amounts of energy to 

c/3 ! large distances from the source, and their interaction with the ambi- 

' ent medium has a crucial effect on the evolution of the system. The 

propagation of the jet is characterized by the formation of a shocked 
"head" at the front of the jet which dissipates the jet's energy and a co- 

\^ [ coon that surrounds the jet and potentially collimates it. We present 

2^ ' here a self consistent, analytic model that follows the evolution of the 

jet and its cocoon, and describes their interaction. We show that the 
critical parameter that determines the properties of the jet-cocoon 

(^ ' system is the dimensionless ratio between the jet's energy density and 

the rest-mass energy density of the ambient medium. This parameter, 
together with the jet's injection angle, also determines whether the 
jet is collimated by the cocoon or not. The model is applicable to 

^ ■ relativistic, unmagnetized, jets on all scales and may be used to deter- 

mine the conditions in AGNs jets as well as in GRBs or microquasars. 
It shows that AGN and microquasar jets are hydrodynamically colli- 
mated due to the interaction with the ambient medium, while GRB 
jets can be collimated only inside a star and become uncollimated once 
they breakout. 
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1 Introduction 



Relativistic jets are observed in many astrophysical systems which host com- 
pact objects, such as radio galaxies, gamma-ray bursts (GRBs) and micro- 
quasars. These jets appear in a variety of lengths, durations and energy 
scales, and propagate in different types of media. However, despite the dif- 
ferences in their characteristics, which are many orders of magnitudes apart, 
the interactions of the jets with their surroundings lead to similar results, 
such as energy injection into the ambient medium and its feedback on the 
jet. Understanding these common phenomena can provide valuable insights 
on puzzles such as what type of medium can collimate the jet and how much 
energy is injected into the surrounding matter. These in turn can be used 
to study the heating of the interstellar and intergalactic medium by active 
galactic nuclei (AGNs) and the fate of jets and stellar envelops of collapsing 
massive stars during gamma-ray bursts. 

The interaction of a jet with the external medium was studied exten- 
sively in many different scales, using analytic and numer ical methods. Ini- 
tially it was studied in the context of radio loud AGNs (JBlandford &: Rees 
1974 IScheuerlll974l ). They showed that the propagation of the jet generates 
a double bow-shock structure at the head of the jet. Energy and matter 
that enter this structure are pushed aside due to a high pressure gradient 
and create a hot cocoon around th e jet. The cocoon, in tur n, applies pres- 
sure on the jet and compresses it. iBegelman &: Cioffil (119891 ) calculated the 
head's velocity and the cocoon expansion rate, assuming that the cocoon and 
the head are supported by the ram-pressure of the ambient medium. Their 
calculations assumes Newtonian head velocities and rely on the knowledge 
of the cross-section of the head, which they estimated phenomenologically 
using the size of the r adio l obes observed around jets of radio loud AGNs. 
Meszaros &: WaxmanI (1200 ll ) calculated the Lore ntz factor of a relativistic 
head assumi ng a conical jet (i.e. , no co llimation) . JMatznerl ( l2003l ) extended 
the model of JBegelman fc Cioffil (11989! ) to include bot h Newtonian and rel- 
ativist ic head velocities, assuming a conic al jet as in JMeszaros fc Waxman 
(1200 ll ). Later, JLazzati fc Begelmanl ( 120051 ) used this model to calculate the 
opening angle and the properties of the jet at breakout from a surface of a 
collapsing star. They accounted for the reduction in the opening angle of the 
jet due to collimation, but ignored the dissipation of energy in shocks that 
form within the jet as a result of such a collimation. This led to an incorrect 
dependence of the opening angle on the cocoon's pressure. An attempt to in- 



elude the effects of these shocks was made by lMorsony. Lazzati.fc Begelman 
(l2007t ). 

Apart from the analytical studies the numerical approach was also given 
a lot of attention over the years . Simulati ons were conducted in the con- 



text of extra galactic jets (e.g . iMarti et al.l Il995l . 119971 : lAloy et al.l Il999 : 



Hughes. Miller. Sz Duncarul2002[) as well as GRB jets that propagate inside a 
star [e.g. IZhang. Wooslev. fc MacFadvenll2003l : iMorsony. Lazzati.fc Begelman 
20071 : iMizuta fc Aloyl l2009l ) . All these works showed the basic features dis- 
cussed above, i.e. the formation of a jet head and a hot cocoon. A coUimation 
shock at the base of the jet is also evident in the simulations, whenever the jet 
is collimated by the cocoon. However, despite the extensive numerical and 
analytical efforts, to this date there is no simple analytic description of the 
time evolution of the jet-cocoon system. As a result, there is no quantitative 
understanding of the collimation process, and its effect on the jet angle. 

The goal of this work is to provide a simple self consistent, analytic de- 
scription, for the time evolution of an unmagnetized relativistic jet and the 
cocoon that forms around it, when it propagates in a density profile that is 
suitable for most astrophysical systems. The key point in our analysis is the 
treatment of the collimation shock that forms at the base of the jet. This 
shock is an inevitable consequence of collimation in a supersonic jet, as it 
dissipates part of the jet's energy, generating the pressure needed to coun- 
terbalance the cocoon's pressure. It is crucial for the proper modeling of the 
system, since it sets the width of the jet and cont rols its propagation veloc - 
ity. We base ou r solution on the earher analysis of iBegelman fc Cioffil (119891 ): 
Matznerl (|2003|). and incorporate the geo r netry of the collimatiori shock as it 
is analyzed by iKomissarov fc Falld (jl997l ) ; iBromberg fc LevinsonI (120091 ) . We 
find that the jet evolution can be of two types: collimated and un-collimated 
according to the strength of the jet - cocoon interaction and the collimation 
shock. We quantify the transition between these regimes and derive models 
for the jet and the cocoon evolution in each regime. We test the validity of 
our approximations by comparing the results from our model with output of 
various numerical simulations. 

The paper is constructed as follow: In §2 we present a general nontech- 
nical description of the system, its ingredients and the different stages in its 
evolution. In §3 we present our model and derive the criterion that distin- 
guishes between the two collimation regimes. We discuss the geometry of 
the collimation shock in each regime and analyze how it affects the temporal 
evolution of the system. We also provide approximated analytic solutions 



in each of the regimes. We compare our analytical model to hydrodynamic 
simulations in §4. In §5 we briefly consider jet collimation in various astro- 
physical environments. Thorough study of specific systems is lengthy and 
will be presented in following papers. In §6 we summarize the main results. 
The full list of results, that specifies the behavior of the physical quantities 
in the entire parameter space is summarized in table 1 and in appendix B. 

2 A schematic description of the system 

Consider an un-magnetized, sufficiently faslo relativistic jet with a luminosity 
Lj, injected at a fixed opening angle 6o, into the surrounding medium. We 
assume that the entire system (jet, cocoon and ambient medium) is axisym- 
metric and use cylindrical coordinates {z,r). We further approximate the 
ambient medium density, pa, to depend only on the hight, i.e. pa{z). This is 
a good approximation even in cases where pa is spherically symmetric (e.g., a 
stellar envelope), since the opening angles of the jet and the cocoon are small 
whenever the the ambient medium plays an important role and collimates 
the jet. 

The jet propagates by pushing the matter in front of it, leading to the 
formation at the jet's front of a forward shock and a reverse shock that are 
separated by a contact discontinuity. We refer to this structure as the jet's 
head. Matter that enters the head through the shocks is heated and flows 
sideways since its pressure is higher than that of the surrounding matter (see 
fig. 1). This leads to the formation of a pressured cocoon around the jet. A 
contact discontinuity divides the cocoon into an inner, light, part containing 
the jet material which crossed the reveres shock, and an outer part of the 
heavier shocked medium. The pressure in both sides of the discontinuity 
is equal, but since the plasma is much more tenuous in the inner cocoon, 
the sound speed there is much larger than in the outer cocoon. This allows 
a causal connection along the cocoon, leading to pressure equilibration in 
the vertical direction and a more uniform distribution of the energy (see 
further discussion in|3]). Since our analysis needs only the cocoon's pressure 
we disregard in the following the cocoon's inner structure. If the cocoon's 
pressure is sufficiently high, it collimates the jet and reduces its opening 
angle. This changes the jet's propagation velocity and the energy flow into 



^Defined later at §3.1, eq. (HI]). 



the cocoon. We classify two regimes of the system: coUimated and un- 
coUimated, according to how strongly the jet is affected by the cocoon's 
pressure. Fig. ([T]) depicts a schematic description of the system in these two 
regimes. 

In the collimated regime (fig. [H right panel), the cocoon's pressure col- 
limates the jet and reduces substantially its opening angle. Since the jet is 
supersonic, the collimation leads to the formation of an oblique shock at the 
base of the jet. This shock deflects the jet flow-lines and generates a pres- 
sure that counterbalances the cocoon's pressure. To maintain the required 
pressure the shock curves toward the jet's axis, until it converges at some 
altitude, above which the collimation is complete (see §3.1). The geometry 
of the collimation shock sets the jet's cross-section at the head to be much 
smaller than the cross section of an un-collimated jet. Consequently the jet 
applies a larger ram-pressure on the head and pushes it to higher velocities. 
The faster motion reduces the rate of energy flow into the cocoon. At the 
same time the cocoon's height increases at a faster rate, resulting in a larger 
volume and a decrease in the cocoon's pressure. There is a limit to the head 
velocity above which the cocoon's pressure is too low to effectively collimate 
the jet. We show that this occurs when Lj/^z^paC") — 6q , corresponding 
to a head's Lorentz factor Th — 9q . Therefore for typical initial opening 
angles ^o > 1°, the head velocities in the collimated regime can be at most 
mildly relativistic. 

The un-collimated regime is characterized by larger values of T^ and a 
cocoon pressure which is insufficient to collimate the jet in an appreciable 
amount. The jet remains conical to a good approximation and the collimation 
shock remains at the edges of the jet and does not converge onto the jet's 
axis. This results in a coaxial jet structure composed of an inner fast spine 
surrounded by a denser layer of the shocked jet material, having a lower 
Lorentz factor (Fig. [T], left panel). When T^ > 9q^, the head moves so 
fast that different parts of the jet's head become causally disconnected and 
energy can flow into the cocoon only from a small region of the head. This 
further reduces the cocoon's pressure. At even larger Lorentz factors the 
reverse shock at the head becomes weak, and it no longer affects the jet, 
which can be considered as propagating in a vacuum. The forward shock 
continues, however, to gather matter from the ambient medium in front of 
the jet and to accelerates it to a Lorentz factor similar to that of the jet. 
A small fraction of this shocked matter continues to stream into the cocoon 



and feeds with relativistic particles. 

Following the above description we divide the system into five main ele- 
ments (see figure [T]): un-shocked jet, shocked jet (by the coUimation shock), 
jet's head, cocoon (containing an inner light part and an outer heavy part) 
and ambient medium. The system dynamics is determined by the following 
relations between these region, (i) The jet's head velocity is set by balancing 
the ram-pressure applied on the forward shock (by the ambient medium) with 
the ram-pressure applied on the reverse shocks (by the shocked or unshocked 
jet, depending on the coUimation regime). The head velocity determines the 
cocoon height and the energy injection rate into it. The jet head has the 
highest pressure in the system, (ii) The pressure in the cocoon is set by its 
size and by the energy injected into it from the head, (iii) The velocity of the 
shock that infiates the cocoon into the ambient medium is set by the balance 
between the ram-pressure of the ambient medium on that shock and the co- 
coon pressure, (iv) The pressure in the shocked jet is equal to the cocoon 
pressure. The coUimation shock structure is set to build up this pressure in 
the jet. This structure, in turn, determines the jet head cross section and 
thus the ram pressure applied on the reverse shock, (v) The un-shocked jet 
properties are determined by the inner engine. 

3 The jet-cocoon model 

Our model contains the five elements discussed above, where the ambient 
medium serves as a fixed background. Given a jet with a luminosity Lj, 
an injection angle, 6q, and a medium density profile, Pa{z), we calculate 
the time dependent quantities: the head velocity, Ph (predominantly in the 
z direction), the cocoon's pressure, Pc, the cocoon expansion velocity, /3c 
(predominantly in the r direction), and the jet's cross-section S^. We use 
the subindices j, h, c, a to designate quantities related to the jet, the jet's 
head, the cocoon and the ambient medium respectively. The distinction 
between the shocked and un-shocked jet is relevant only in the coUimated 
regime, and we use it when we discuss this regime. The subindex j is used to 
describe general properties of the jet (like the jet dimensions) and when it is 
unimportant which of the two region (shocked or unshocked jet) is considered, 
e.g., the jet luminosity is equal in the two regions and is denoted Lj. 

We begin by discussing the main assumptions of the model. We take Lj, 
9q and p{z) as given, and derive equations for /3ft, P^ Pc and Sj. The later 
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Figure 1: A schematic description of the jet geometry in the two colhmation 
regimes. Left: A coUimated jet; Center: An un-colhmated jet. In both 
panels the basic ingredients of the model are evident: the jet (divided into a 
shocked and and un-shocked part), the jet's head, the cocoon and the ambient 
medium. Also shown are the collimation shock and the contact discontinuity. 
The collimation shock splits the jet to an un-shocked region and a shocked 
region. The contact discontinuity separates the jet material that enters the 
head from the ambient medium. This discontinuity extends to the cocoon 
and divides it to an inner and an outer part. The cocoon expands into the 
ambient medium behind a shock that extends the forward shock at the head. 
All shocks are marked with dashed lines, and the contact discontinuities with 
solid lines. Right: A closeup of the jet's head and the contact discontinuity. 
The matter flows into the head through a forward and a reverse shock, and 
from there to the cocoon as illustrated by the four arrows. 



differs between the two coUimation regimes, and therefore we present the 
relevant solution separately in each regime. 

The head velocity is set by balancing the ram-pressure applied on the 



forward and reverse shocks (JBegelman fc Cioffilll989 



Matzner 



20031): 



p.h^T'^Ti (/3, - P^,y + P,. = PahaTiPf, + P„ 



(1) 



where p, P and h = 1 + AP/pc^ are the mass density, pressure and dimen- 
sionless specific enthalpy of each fluid, measured in the fluid's rest frame, and 
rp is the proper fluid velocity in the frame of the ambient medium. When 
the temperature of the ambient medium is non-relativistic and as long as the 
reverse shock is stron g, Pg and Pj c an be ignored. Under these conditions 
the head's velocity is (JMatznerl 120031 ): 



Ph 



/3., 



1 + L-V2- 



where the dimensionless parameter 

L 



p,h^] 



L, 



Pa 



T^jpaC^ 



(2) 



(3) 



represents the ratio between the energy density of the jet {Lj/T,jc) and the 
rest-mass energy density of the surrounding medium at the location of the 
head. Note that by defining L as such, we implicitly assume that the energy 
lost in the jet due to work against the pressure of the cocoon is negligible 
and that radiation losses are negligible as well. We discuss this assumption 
at length in Appendix A, and show that it is valid as long as the jet injection 
angle, 6o, is not too large. As we show in the following sections, L is the crit- 
ical parameter that determines the evolution of the jet and the cocoon, while 
the combination L9q distinguishes between collimated and un-collimated 
jets. We stress that L may vary with the propagation of the jet, even if the 
jet luminosity is constant, depending on the density profile of the ambient 
medium and on the behavior of the jet's cross section. Thus, a jet can switch 
between the different coUimation regimes, from collimated to un-collimated 
and vice versa. Specifically a collimated jet that encounters a sharp density 

decline may becomes uncollimated. 

In the limits L ^ 1 or L ^ 1, eq. ([2]) can be linearized (JMatznerl l2003l ): 
When L < 1 



/3h ^ L"\ 



(4a) 



corresponding to a non-relativistic head velocity. If 1 ^ L ^ AT^ 

\l"\ (4b) 

which corresponds to a rela tivistic proper velocity of the head, ThPh > 1 (e.g. 



Meszaros &: Waxmanll200ll ). When L > 41^, it follows from eq. ([2]) that 
F/i ~ Tj, implying a non-relativistic reverse shock and a negligible energy 
flow of the shocked jet material into the cocoon. This is equivalent to the 
requirement that the volumetric-enthalpy ratio of the jet and the ambient 



medi um, -^^ > F^, and is similar to the criterion given by ( iSarifc Piran 



paha 



19951 ) for the case of a cold jet. 

The pressure in the cocoon is sustained by a continuous flow of energy 
from the head. At any given moment the total energy in the cocoon is Ec = 
r]Lj{t — Zh/c), where the second term in the brackets represents the fraction 
of the energ y which is carried by the r elativistic jet and hasn't reached the 



cocoon yet (JLazzati fc Begelmarul2005l ). The parameter, 77, varies between 



and 1. It stands for the fraction of the energy that flows into the head and 
enters later into the cocoon. We assume that all the energy in the region of 
the head that is in causal connection with the cocoon flows into the cocoon, 
so 7] is the fraction of that region out of the total head volume. When T^^ 1 
most of the energy that flows out is initially tapped to the bulk motion of 
the fluid and it can't contribute to the pressure. However, later, as the 
matter enters the cocoon it spreads sideways and decelerates exponentially 
to r p^ 1, in a similar manner to a ra pidly spreading jet in a GRB afterglow 
(e.g. iGranot et al.ll2000l : |Piranll2000l ). Therefore practically all the energy 



that flows into the cocoon is available to generate pressure. 

The energy in the cocoon is shared between two parts which are sepa- 
rated by a contact discontinuity: an outer part made of the shocked ambient 
medium with a typical density ~ p^, and an inner part that consists of the 
lighter jet material that crossed the reverse shock and enters from the head 
(see fig. 1). The outer cocoon is supported from the sides by the ram pres- 
sure of the surrounding medium, which balances its pressure. This results in 
a lateral expansion velocity of /3c = \/Pc/ PaC^ which is just below the sound 
speed in the outer cocoon as long as /3c is sub relativistic. In the inner co- 
coon the sound speed is relativistic and typically much faster than /3c. Thus, 
matter on both sides of the contact discontinuity is in causal contact in the 
lateral direction and any pressure difference in this direction is smoothed 

9 



out. In the longitude direction the inner cocoon connects the upper parts of 
the cocoon with lower parts which otherwise would have been out of causal 
contact. This suppresses the pressure drop in the vertical direction, and in 
cases where the head is sub relativist ic leads to a uniform distribution of 
energy in the cocoon. The size of the region connected by the inner cocoon 
is determined by the survival length of the contact discontinuity, which is 
set by the growth rate of Kelvin-Helmholtz instabilities. The small shear- 
velocity differences on both sides of the contact render the growth rate small 
and stabilize the contact discontinuity. Indeed numerical simulations show 
that the contact discontinuity remains intact for a considerable fraction of 
the cocoon's length, and that the pressure drop is suppressed (see §4). We 
stress that the uniform pressure approximation should also hold rather well 
in cases where Pc develops some vertical gradient since the conditions in the 
jet's head depends only weakly on the pressure gradient. 

To calculate the cocoon's pressure we make the following approximations: 
First we take the cocoon's shape to be a cylinder with a hight Zh = J Phdt and 
a radius re = J f3cdt. Second, based on the discussion above, we approximate 
the energy density to be uniformly distributed within the cocoon. Third, 
over a wide range of parameters and especially in the systems relevant for 
us, the pressure in the cocoon is radiation dominated, thus we use in our 
analysis an adiabatic index of 4/3. Under these approximations the cocoon's 
pressure satisfies: 

To calculate the transverse velocity of the cocoon we take the average density 
of the medium pa = 1/Vc{z) f pa{z)dV, where Vc{z) is the volume of the 



cocoo n, and obtain the lateral expansion velocity: (e.g. iBegelman &: Cioffi 



19891): 



Pa{Zh)c^ 



^c = xI-t:^:^. (6) 



Approximating / I3hdt ~ (3ht and / f3cdt ~ /3ct and substitute into eq. (J5]) we 
obtain the cocoon's pressure: 

Pc = s,[^]'^'L-V4r\ (7) 
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where S^ is a constant of order unity that depends on the density profile 
of the ambient medium, and on the colhmation regime. The value of S^ in 
each regime is given in the appendix, where we list the full time dependent 
solutions of eqn. ([2]), (jS]) and (jH]) in the limits: L -C 1 and L ^ 1. 

3.1 The coUimated regime {L < 6q ). 

In the collimated regime the pressure in the cocoon is sufficiently strong to 
compress the jet and form a strong oblique shock at the base of the jet. The 
jet's head is then at full causal contact and therefore rj = 1. The collimation 
shock divides the jet into a pre-shocked and a shocked region. We designate 
them by subindices jO and jl respectively (See fig. [2]). The fiow-lines in area 
jO are radial and they encounter the shock at an angle "^{z), defined as the 
angle between the original direction of the fiow-line and the shock surface. As 
the fiow-lines pass through the shock, they are defiected toward the jet's axis. 
Some of the fiow energy is dissipated generating the pressure that supports 
the jet against the compression of the cocoon. The shock's geometry is set 
by the balance between the cocoon's pressure and the u pstream momentum 



fiux, normal to the shock, that en ters at the upstream ( iKomissarov fc Falle 



19971 : iBromberg fc LevinsonI 120071 ): 



Pjoc'hjoT'^o(3-o sin' ^ + P,o = Pc- (8) 

Since the jet is relativistic and its ram pressure is much larger than its internal 
pressure, Pjo can be neglected on the LHS of eq. ([8]). The radial geometry 
of the unshocked fiow-lines implies that pjoc^hjoT^Q/3'jQ oc z~^. Therefore 
to keep a uniform pressure at the downstream \E^ must increases with z. 
Consequently, the shock curves toward the axis until it converges on it at 
some point. In the small angle approximation, to a first order, sin\E^ = 
(y — ^); where r^ is the cylindrical radius of the collimation shock (see 
fig. 2). In this limit eq. (Ml becomes a first order OD E whose solution is 



(IKomissarov &: Fallelll997 : 



Bromberg Sz LevinsonI l2009l ): 



9oz (1 + Az,) - A9oz^, (9) 



where A = \ j^§^-, and z^^ is the hight where the jet is first affected by the 

cocoon and the collimation shock forms. As long as the pressure in the jet 
is larger than the pressure in the cocoon, the jet is unaware of the cocoon. 
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As the jet accelerates pressure is gradually converted to kinetic energy until 
eventually it reaches a point where Pjo{z) = Pc- At this point the jet's 
compression by the cocoon becomes significant. Therefore if the jet-medium 
interaction begins at the injection point z^^ can be tracked to the point where 
PjQ = Pc. Alternatively, if the jet is injected into a cavity of radius i?*, 
Zj, = Taax{R^, z{Pjo = Pc)}- Note that once the shock forms, its geometry 
depends only on the jet luminosity and it is the same for an accelerating or 
a non-accelerating jet. The shock converges at r^ = 0, namely at 

z = A-^ + z,. (10) 

As long as z^ < A~^ it can be neglected from eqn. flQfTOl) . which implies 
that i?* must be smaller than A~^ as well. When 2;* > R^, we can use the 
jet luminosity at 2* which holds Lj > in PjoT'jqI3joz19qC together with the 
condition that Pjq = Pc and extract z^. Substituting A, we then get that 
z^ < A-^ if 

Tjoiz*)Pjoiz*) > Oo'- (11) 

We define a "sufficiently fast" jet as one that satisfies this condition. In 
this situation we can ignore 2* as long as i?* < A~^. Since we work in the 
limit where ^0 is small, eq. (11) also implies that the jet is relativistic and 
therefore in the following we approximate 13 jq = 1 and /3ji = 1. Note that as 
long as Tjq{z) < 9^^ the jet expands sideways under its own pressure and its 
opening angle increases as it propagates. In such a case it is meaningless to 
discuss a constant initial opening angle, and our model does not hold. The 
jet stop expanding when the angle of the fiowlines is of the order of F,^^, one 
over the injected Lorentz factor. Beyond this point it corresponds with our 
condition for neglecting z*. In this case we can approximate 6*0 to be ~ Tq^ 
and neglect 2;*. 

To estimate the jet's cross-section we approximate the jet to be conical 
up to the point where collimation has a sizable effect on the jet's geometry, 
which is where the collimation shock is roughly parallel to the z axis. From 
eq. iQ this occurs at z/2. Above this point the jet is collimated and since 
the cocoon pressure is roughly uniform, the jet cross-section doesn't change 
much, and it can be taken to be constant. In this approximation the cylindric 
radius of the jet is rj{z > ^) ~ z9q/2 and its cross-section is: 

EAz > 5) ~ -TTzX - ^- (12) 

^^ ^40 4cPc ^ ' 

12 



The luminosity above z is 



h ^ 4P,r^iS,c. 



(13) 



Substituting Lj in eq. ( !T2|) we obtain the Lorentz factor of the shocked jet 
material above z: 

r,i(^>i) = l (14) 

This implies that, as long as the jet is sufficiently fast at the injection point, 
the Lorentz factor of the collimated jet depends only on the injection angle, 
and it is independent of the initial Lorentz factor. 
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Figure 2: A schematic description of the collimation shock in the collimated 
regime. The collimation shock is marked in a dashed black line, and it 
separates the jet into an un-shocked region (jO) and a shocked region (jl). 
The jet flow-lines, marked with light blue arrows, are initially radial. They 
intersect the shock at point (r^,^) and form an angle '^ with the shock's 
surface. Downstream of the shock the flow-lines are deflected toward the 
jet's axis. The shock crossing point of a speciflc flow-line is illustrated. The 
shock converges onto the jet's axis at z. Above z the jet maintains a constant 
cylindric radius, r^. 

The remaining system parameters are calculated by substituting eq. flT2|) 



13 



for the jet's cross section into equations (3-7). Eqn. (|3]) and (j7]) yields: 

f - £_^ ~ A+-4/5r 2/5 2/5^-8/5 .-. .s 

and 

P, ^ t-^/^Lf Va^'^o^'- (16) 

The head's position is calculated using Zh ~ Pht-, and the linearized expression 
for Ph given in eqn. (4a,b): 

L < 1 

Zh-^ { (17) 




Taking the cocoon's cylindric radius r^ ~ t\J P^Pa and the jet's cylindric 
radius r^ ~ Oq^JLj/cPc^ we calculate the opening angle of the cocoon and 
the jet respectively: 
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and 

f ZV4^2 ^ ^ << 1 

(19) 

Note that if the head is non-relativistic, the cocoon's aspect ratio is constant 
and it equals the jet's injection angle, up to a constant of order unity (see 
Appendix B). 

We can now infer the parameter regime for which the jet is collimated 
by the cocoon. The jet is considered collimated if the collimation shock 
converges below the jet's head. Therefore, the transition to the un-collimated 
regime occurs when z = Zh- We find that this equality can take place only 
when L > 1, thus we substitute t = Zh/c in eq. (fT6l) and get that for Zh > z, 
Pc > PaC^^Q . Substituting Pc in eq. ( !T5|) . we find, up to a constant of order 
unity, that the condition for collimation is: 

L < e-"\ (20) 
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The numerical factor missing in this equation depends on the density profile 
(see appendix B). This condition corresponds to a limit on the head's Lorentz 
factor 

r. < e,'/\ (21) 

—1/3 

If the head accelerates beyond 6q the collimation shock fails to converge 
below the head, and the jet can be approximated as having a conical shape. 
The jet can shift from a coUimated to an un-coUimated state and vise 
versa, depending on the density profile of the ambient medium. According 
to eqn. f lTSllTTl) . L oc Pjpa oc {z-^p-^f, where 5 = 1/3 if L < 1 and 5 = 1/5 
if 1 ^ Z < 6'q' . If the density profile is steeper than z~'^, then L increases 
with z and the jet's head accelerates. In such a case 6j = Oo^ increases 

with time. When L ~ ^(^ , the collimation shock converges at the head, 
and 9j = 6*0/2. At higher values of L the shock opens up and the jet becomes 
conical. In the limit of pa oc 2;"^, L is constant and the jet's head velocity is 
constant. Eq. (1191) shows that in this limit 6j is constant as well, implying 
that Tj oc Zh- This makes an interesting case where the jet expands sideways 
at the same rate it propagates upward. Although the opening angle of the 
jet, 6j <C ^0; it remains constant like in the case of a conical jet. 

3.2 The un-collimated regime 

When L ^ 9q {Th ^ ^q ) the pressure in the cocoon is too weak to 
significantly affect the geometry of the jet. In this case the jet remains 
conical to a good approximation and 

S,(^,) = nzlel (22) 

The collimation shock remains at the edge of the jet, resulting in a coaxial jet 
structure of a cold, fast, inner spine surrounded by a hotter and denser layer 
of the shocked jet material moving with a lower Lorentz factor (see fig. 1). 
The layer of shocked material becomes thinner and thinner at higher values 
of L. In this regime most of the jet's plasma streams freely all the way to 
the head (region jo extends to the head), and dissipates all its energy in the 
reverse shock. The jet is therefore cold below the head with some brightening 
at its limb, as opposed to a coUimated jet which is hot below the head, since 
its plasma is first shocked much closer to the base by the collimation shock. 
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Substituting Sj in eq. ([3]) and taking Zh = ct gives: 



^3 



TiOlp^t^C^- 



(23) 



This implies that also in the un-coUimated regime pa ~ z~'^ remains the lim- 
iting profile which distinguishes between an accelerating and a decelerating 
head. 

The cocoon has some distinctive properties which can be classified to 
three subcases, according to the values of L and r/, the fraction of energy 
that fiows from the head into the cocoon (see also table 1). In each subcase 
the conditions in the head are different, and this affects the amount of energy 
that fiows out of the head: 

1) e^^'"^ < L < ^^^ {e^^^^ < T/, < ^^1): The head is sufficiently slow. 
Different regions of the head are in a causal contact with each other 
and all the energy in the head can flow into the cocoon, thus 77 ^ 1. 
The cocoon's pressure is calculated by substituting L from eq. ( 123|) in 
eq. ([7]) and using Zh = ct: 

P, ^ t-V^Lf p3/4^1/2^3/4 ^ Z^/^^oPaC^ (24) 

Under this pressure the temperature at the outer cocoon remains sub- 
relativistic and the different parts of the cocoon maintain causal con- 
nection in the lateral direction. As L increases Pd Pa(? grows and it 
approaches unity as L — )■ Q'^'^ . In this limit the pressure becomes mildly 
relativistic and it pushes the edge of the cocoon to a velocity /3c —^ 1, 
which is above the local sounds speed, c/v3. This results in a loss of 
causality in the transverse direction which implies that the approxima- 
tion of a uniform pressure no longer holds. In addition the cocoon's 
aspect ratio approaches unity and it can no longer be considered as 
cylindric. Therefore above the limit of L ~ B'^'^ our model can only 
provide the total energy in the cocoon. This has no effect on the con- 
ditions in the jet and the jet's head, which are well described by our 
model, since in the un-collimated regime the jet is not sensitive to the 
pressure in the cocoon. 

2) 6'o''^ < Z < r|o (6'o'^ <^h < Tjo): At this limit the reverse shock 
is still strong and pressure in the head is large, but the jet's head loses 
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causal connectivity in the transverse direction. Energy can flow to 
the cocoon only from a thin annulus on the edge of the head with an 
opening angle F^^, which corresponds to a fraction r] = 2{Th9o)"^ of 
the total energy that enters the head. Most of the energy that flows 
into the head remains trapped and accumulates at a rate of LjTj^^. The 
total energy that flows into the cocoon can thus be estimated as 

E, = vLjTft ~ L'/%Pach\ (25) 

where we substitute Lj from eq. (l23l) and used the relation L ~ 4r^ 
from eq. (4b). Since (3^ = 1, the cocoon has a spherical rather than 
cylindrical shape and it occupies a volume Vc — (ct)^. 

3) TjQ ^ L [Th — Tjo): In this regime the reverse shock becomes 
Newtonian while the forward shock remains relativistic and moves with 
a Lorentz factor ~ Tjq. This leads to different sound speeds in the head, 
above and below the contact discontinuity. Below, within the shocked 
jet material, the sound speed -C c. Above, within the shocked medium, 
the sound speed is c/-\/3. As a result the energy that flows into the 
cocoon comes mostly from the ambient medium part of the head and it 
flows into the outer cocoon, while the inner cocoon which is fed by the 
shocked jet material in the head becomes insignificant. To measure the 
energy that enters the cocoon we first estimate the energy per unit time 
that flows into the head through the forward shock: Eha = Pa^'jo^o'^'^c^- 
The total energy in the cocoon is therefore: 

E, = rjEhat - TjoeoPat^c^ (26) 

where r] = 2(Tjo9o)~^ in this case. 

Table. [1] below summarizes the different characteristics of the jet-cocoon 
in the four different collimation regimes. 

3.3 Comparison with previous analytical works 

The propagation and interaction of a jet with an ambient medium was studied 
in the past in various parameter regimes. Here we briefly comment on some 
relevant works and discuss their compatibility to our results. 

The main new feature introduced in this work is a proper closure of the 
set of equations using the conditions at the collimation shock. This allows 
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us to determine if the jet is coUimated or not, given its luminosity, initial 
opening angle and the ambient density. When the jet is coUimated it also 
allows us to determine its cross section, Sj. Earlier works differ from ours 
in the way they deal with this issue. Most ignore this coUimation shock and 
assume a specific value or functional form for one of the variables (e.g. the 
jet cross section), obtaining a closure in this way. Such solutions typically 
agree with ours if we substitute the value we obtain from the full system of 
eq uations instead of the free variable in these solutions. 



Begelman fc Cioffil ( 1l989l ) analyzed the propagation of a non-relativistic 



galactic jet in the inter galactic medium, assuming that the jet is coUimated 
and that it has a constant velocity Vj < c. When using our expression for Sj 
(eq. [T2|) in their solution ( and taking the li mit where the Vj — ?■ c), their results 
agree wit h ours for L -C 1 . Matznen ( 20031 ) extended this model to relativistic 
jets. Like iBegelman &: Cioffil ( 19891 ) he did not model the collimation of the 
jet and used 6j (which is related to Sj) as a given parameter. This solution 
is consist ent with ours (f or L ^ 1) upon substitution of 9j from our solution. 
However, iMatznerl ( 20031 ) assumed that the jet is in fact conical, with 9j = 9q, 
which is inconsistent in this limit of L <C 1. 



Lazzati fc Begelman! (120051 ) . considered the collimation of the jet by the 



cocoon's pressure. However they ignored the dissipation in the collimation 
shock. Instead they assumed that the jet material expands adiabatically, 
leading to a relation Tj oc S- . As we show in sec. 3.1 (eqn. [T2lfT4|) . 
the collimation shock renders this relation invalid. Consequently, the so- 
lution they obtained differs from ours in all regimes. Their solution has 
a smaller opening angle and a large r value of F,- at breakout. Later on. 



Morsony. Lazzati. &: BegelmanI ( 120071 ) attempted to calculate the geometry 
of the collimation shocks. But they have used an incorrect expression for 
the momentum flux that crosses this shock, which led to a shock that never 
converge to the axis (z — »■ oo). 



Finally, iMeszaros fc Waxmanl ( 1200 ll ) analyzed the propagation of an un- 
coUimated GRB jet in the outer envelope of a red supergiant. They consid- 
ered only the properties of the jet's head, ignoring the surrounding cocoon. 



Their solution for T^ is valid for ^n <^ L <^ Ff 



Table 1: The system's characteristics in the colhmated and un-colhmated 
regimes 
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All quantities are missing order of unity constants of integration over the density profile. In case of 

a power-law density profiles these constants can be calculated analytically and they are given in the 

appendix. 

f Tj is the jet Lorentz factor just below the head. 

'f'l' This quantity represents the total energy in the cocoon divided by the volume Vc = (ct)^ 

and not the real pressure, which is not calculated in this regime. 

4 Comparison with numerical simulations 



Jet simulations have been carried out extensively by various authors. We turn 
now to compare our results with two recent numer ical simulations. First, we 
consider the simulation by iMizuta fc Aloyl ( 120091 ) who modeled the propa- 
gation of a relativistic jet in the envelope of a massive st ar. To model the 
star th ey used a numerically calcu lated density profile fr om lWoosley fc Heger 
(l2006l ) (model number HE16N in lMizuta fc Alovll2009h . This star has a ra- 
dius of ~ 6 X 10^*^ cm, and it has a density profile which can be divided into 
three parts. The inner part (up to ~ 1.2 x 10^*^ cm) has an average powerlaw 
profile with an index a = —2.5. Above it (up to 4 x 10^° cm) the profile is 
steeper with an averaged index a ^ —4.5, and it drops sharply from there 
to the edge of the star. The jet in the simulation is hot and it is injected 
with a Lorentz factor, Fq = 5 into a cone of an opening angle Oq = 5° and 
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an initial altitude of 2* = 10^ cm. Since 9q < Tq^ the jet goes through an 
initial transitory phase where it's opening angle increases until it stabilizes 
when Ojo{z) ~ TJq{z) ~ 10°. We therefore take in our calculations Oq = 10°. 
Figure (I3l) shows two snapshots of the jet and the cocoon at times 0.8 sec 
and 1.2 sec after the injectioiij. The times were chosen so that the jet is 
far enough from the injection point and the opening angle at z* is already 
stable. The color codes represent equal values of the pressure normalized by 
1/c^ (left side) and density (right side). Our calculation of the jet, the cocoon 
and the collimation shock is drawn in black lines on each panel. Due to the 
relatively large 2;* (see below), we did not use the analytic formulae, instead 
we have integrated eqn. (2,5,6,10,12) numerically to obtain a solution. It can 
be seen that in a density profile of a = —2.5, the approximation of a uniform 
pressure in the cocoon fits well with the pressure profile in the simulation. 
Our results agree with the numerical results to a good accuracy. 

Figure (jlj) presents two snapshots of the jet and the cocoon in the region 
with the steeper density profile, just before the sharp drop^. The steeper 
gradient in the stellar density profile leads to an acceleration of the jet's head, 
which changes the rate of energy flow into the cocoon. Nevertheless it can 
be seen that the pressure remains uniform to a good approximation through 
most of the cocoon's hight, due to the survival of the contact discontinuity 
that keeps the inner cocoon intact. The pressure calculated by our model fits 
well with the simulated pressure in the upper, uniform, part of the cocoon. 
At ~ 1.2 X lO^'' cm the contact discontinuity becomes unstable and below this 
point the light matter from the inner cocoon mixes with the heavy material 
in the outer cocoon. As a result the matter in the bottom part of the cocoon 
is no longer in contact with the matter in the upper pars and a pressure 
gradient is formed. As a consequence the collimation shock converges at a 
lower altitude than what we calculate (at ~ 4 x 10^ cm), which leads to a jet 
that is narrower at its base. But due to the negative pressure gradient the 
jet widens to a similar size as calculated by our model, and the jet's head 
has similar properties (width and velocity) as in our calculation. The region 
where the contact is destroyed corresponds to a sharp step in the stellar 
density profile that is located at ~ 1.2 x 10^° and separates the shallow 
profile of a = —2.5 from the steeper one with a = —4.5. At this step the 
growth rate of instabilities on the contact discontinuity increases. 



^'^ The snapshots are a curtesy of A. Mizuta, and are taken from a simulation published 
Mizutafc Alwl (|2009t ). 
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Figure 3: The calculated size of the jet, the coco on and the coUimation shock 
drawn in black lines on top of a simulated jet by lMizuta fc Aloyl (120091 ) . The 
left (right) panel shows a snapshot of the jet and the cocoon after 0.8(1.2) 
sec. Color codes are of equal pressure divided by c^ on the left, and equal 
density on the right. The black arrow shows the average value of the cocoon's 
pressure from our calculation. Figures are at curtesy of A. Mizuta. 



We also compared our model with the simulations of lZhang. Woosley. fc MacFadyen 



(120031 ). They examined the consequences of changing the injection angle and 
the luminosity of the jet on the propagation in a stellar mantle. Table [2] 
shows the initial parameters and the breakout times of the jet from the stel- 
lar surface in the three cases examined. In the right column we added the 
breakout times calculated by our analytical model, using the same stellar 
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Figure 4: The calculated size of the jet, the coco on and the collimation shock 
drawn in black lines on top of a simulated jet by lMizuta fc Aloyl (20091). The 
left (right) panel shows a snapshot of the jet and the cocoon after 1.6(2.4) 
sec. See fig. Elfor further details. Figures are at curtesy of A. Mizuta. 



density profile and initial jet parameters . It can be seen that our results 
agree to within 10% with the simulated ones. 

When comparing our results with numerical simulations we should recall 
that the these simulations usually use a relatively large value of z^, the ra- 
dius in which the jet is injected to the stellar envelope. Such a jet is initially 
wider than a similar jet that is injected at a smaller radius and at this stage 
it propagates slower. Eventually the solution converges but the resulting 



breakout time is longer by 






which becomes significant for 
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z* ^ 10"^-R, where R is the stellar radius. The values of z^. in jet simula- 
tions are usually larger than that and therefore we integrate numerically eqn. 
(2,5,6,10,12) to calculate the breakout times of the jet presented in table 2. 



Table 2: Breakout Times 



Model" Lj/lQ^^erg Oq t' [sec] t^ [sec 



JA 
JB 
JC 



1.0 
1.0 
0.3 



20° 

5° 
10° 



6.9 
3.5 

5.7 



6.3 

3.8 

5.2 



° Data from lZhang. Wooslev. fc MacFadvepI (i2003l) 

^ The breakout time according to the simulation 

'^ The breakout time according to the analytic calculation. 



Though our analytical model greatly simplifies the conditions in the jet 
and the cocoon, it shows a remarkable agreement with the results of the sim- 
ulations presented above, which is slightly better than what we would expect 
in the general case. Generally we expect an order of unity agreement for all 
density profiles that are not too steep, so that the forward and reverse shock 
of the head remain in causal contact and that the approxim ation of uniform 
cocoo n pressure is reasonable. As we see in the comparison to lMizuta fc Aloy 
(120091 ) even a power-law density gradient as steep as pa oc z~'^'^ satisfies these 
conditions. These conditions break down however in a very steep density pro- 
file, such as the one at the edge of a stellar envelope, where the forward shock 
accelerates and lose causal contact with the jet. 



5 Jet collimation in astrophysical environments 



4/3 



According to our model when L > 9^ the jet is uncollimated, and it has 
a conical shape with 9j = Oq. It implies that in this regime L 



-4/3 



■KZ^elpaC^ ' 



When L < 0Q , on the other hand, the jet is coUimated and 9j < 9q. We 



can therefore formulate the condition for collimation as: 



T^ZlpaC 



3<^0 



2/3 



(27) 



and evaluate the conditions in the jet that lead to its collimation in different 
astrophysical media. 
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The conventional view ab out long GRBs is that they originate from mas- 
sive stars that collapse (e.g. IWoosleyl Il993l ). The jet propagation in this 
model is characterized by two stages: pre-breakout and post-breakout. Prior 
to the breakout the jet propagates inside the star which is considered to be 
massive (M ^ lOMp)) a nd rather comp act {R ~ Rq) Wolf Rayet star (e.g. 



Woosley fc Hegerl l2006t ICrowtherl 120071 ) . The density profile in the stellar 
envelope can be approx imated as pa = p{z/R)~"', with 2 < a < 3 (e.g. 
Matzner fc McKeelll999l ). and p = ^^MR~^ is the average density. Substi- 
tuting this in eq. 027p we get that the jet is collimated if: 



^'<^-''&'^"{iTm 



M 



2/3 / 

VIOMq 



erg/s. {21 



The luminosity we measure in GRBs reflects the luminosity of the jet after it 
breaks out. We take here the simple approach that this luminosity is constant 
over time. Under this assumption we can use the observed jet luminosities 
to estimate the collimation of the jet inside the star. Correcting to a typical 
beaming angle of 5° — 10°, the observed values give Lj < 10^^ ergs/s. In 
addition, as we show below, the measured values of the beaming angle of 
the GRB are similar to the size of the injection angle at the base of the 
jet, Oq. Substituting this in eq. fl28|) we get that GRB jets are collimated 
before they breakout. The propagation of the jet inside a stellar envelope, 
an d the implications on G RB observations are discussed to a greater extent 



m 



Bromberg et al.l (120 111 ) 



Once the jet breaks out its opening angle, 6j, can be measured, for exam- 
ple by identifying a "jet break" in the afterglow lightcurve. Therefore it is 
better to express the condition for collimation (eq. [27|) in terms of 6j rather 
than ^0- Generally 6j ^ 6q and 6j = 6q if the jet is uncoUimated. Therefore 
if L > ^7^/^ it follows that L > 9'^'^. On the other hand if L < 9'^'^ then 

necessarily L < 9q (otherwise the jet is uncoUimated and 9j = 9q). Thus 
eq. fl27j) can be used with 9j replacing 6*0. 

The medium outside the star can either be a dense wind ejected from 
the surface of the star, having a typical density profile pa = a*/z'^, where 
a* ~ 5 ■ 10^^ g/cm, or it can be the constant density ISM. If the jet breaks 
out into a stellar wind environment it is collimated if: 



Lj < 10 



43 



10^ 



2/3 



5 ■ 10^^ g/cm 



erg/s. 



(29) 
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If alternatively it propagates in the ISM it is coUimated if 

Since the observed luminosity is orders of magnitude higher than both of 
these limits, the jet becomes un-coUimated when it propagates in these media. 
Therefore, as the jet, which was collimated inside the star, breaks out, it 
rapidly expands sideways and accelerates. Without any external interference 
the opening angle stabilizes when 6j ~ r~j^'^, where Vji is the Lorentz factor of 
the jet at breakout, and since Vji ^ 6'^^ (eq. [T^ we get that the opening angle 
after breakout should be ~ 9q. However the breakout of the cocoon, which 
occurs simultaneously with the jet, limits the sideways expansion. Thus the 
final opening angle of the jet at early times may be smaller than 9q. Note that 
the smaller opening angle only tightens the constraints for collimation given 
in eqn. (1291) and (1301) . At late times, when the cocoon clears out, the jet is no 
longer confined and maintains its initial injection angle. Therefore measuring 
9j at late times gives information about the conditions at the injection site 
of the jet, even before the breakout when the jet was still buried in the star 
and was in fact collimated. 

Microquasars (MC) and X-ray binaries jets, have luminosities of ~ 10^^ 
erg/s, and their opening angles < 10°. Using eq. fl5U]) we get that these jets 
may begin uncoUimated, but beyond ~ 2-10"^ pc they become collimated due 
to the interaction with t he ISM. The collimation of MC and X-ray binaries 



jets was also proposed by iMiller- Jones. Fender fc Nakarl (120061 ). who showed 



that if the values of Tj are low (~ 10) it implies that the jets are collimated. 
Our calculation therefore supports this claim regardless of the Lorentz factor, 
which can not be measured reliably at the present time. 

Powerful Quasar jets extend to distances of hundreds of kpc from their 
sources, and they propagate in a much thinner medium. Their typical power 
is < 10^^ erg/s, and their opening angle < 10°. The density profile in the 
galactic hallo at such distances is usually considered to be isothermal with 
a < 2, and a mass densit y of the order of y 10~^^ — 10~^^ g/cm'^ at a 
distance of ~ 10 kpc (e.g. iBulbul et al.l I2OIOI : ICapelo. Natarajan.fc Coppi 



20101). Substituting these values in eq. (130|1 we get that beyond this distance 



quasars jets are hydrodynamically collimated by their cocoons. 
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6 Summary 

In this work we present an analytical study of the propagation of a relativistic 
hydrodynamic jet in an external medium with a general density profile. The 
interaction of the jet with the medium results in the formation of a shocked 
"head" at the front of the jet, and an over pressured cocoon with a rather 
uniform distribution of energy surrounding the jet. When the pressure of 
the cocoon, P^, is larger than the internal pressure of the jet, it compresses 
the jet and leads to the formation of an oblique shock at the base of the 
jet. Under some conditions, determined in this work, the shock converges 
to the jet's axis, and the jet becomes coUimated. This changes the jet's 
properties, and affects the conditions in the cocoon. Our model follows the 
evolution of the jet, the jet's head, the cocoon and the coUimation shock, 
given three initial parameters: the jet's luminosity, Lj, the jet's opening angle 
at the injection point, ^o, and the density of the ambient medium, Pa{z). We 
determine the condition for coUimation, and provide a self-consistent, time 
dependent, solution to the system's parameters that depends on these three 
initial parameters alone. The main results of our analysis can be summarized 
as: 

• The jet's evolution can be classified into two regimes: a coUimated and 
an un-coUimated, according to the strength of interaction of the jet 
with the ambient medium. The two regimes are distinguished by the 
parameter L and by ^o (or 6*^), which represents the ratio of the jet's 
energy density to the rest-mass energy density of the ambient medium 
at the location of the head. 

• When L < 9q the interaction is strong and the jet is coUimated. In 
this regime the coUimation shock converges on the jet's axis at some 
point, z, below the head. Above z the jet is cylindrical to a good 
approximation and its width is estimated as z6q/2. This width implies 
that the Lorentz factor of the coUimated ejecta satiafies r^i = 9q^. 

• If L -C 1, the head of the coUimated jet is non-relativist ic {/3h -C 1). In 
this limit the cocoon's expansion velocity is proportional to the head's 
velocity, and it has a constant opening angle dc — Oq. 

• When L ^ 6q the coUimation shock fails to converge and it remains 
at the edge of the jet. The jet in this regime in un-coUimated and it 
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remains conical to a good approximation. It has a coaxial structure of a 
fast unshocked spine surrounded by a denser layer of shocked material 
with lower Lorentz factor. 

• When L ^ 6^'^^ V^ > Oq^, and the head is not causal in the transverse 
direction. Here only a fraction ~ (ThOo)^^ of the jet energy flows into 
the cocoon. This energy is enough to produce a mildly relativistic 
pressure in the cocoon, and accelerate cocoon's edge to a velocity (3c ~ 
1. 

• When L > 4TJ^ the reverse shock becomes Newtonian and the head 
no longer affects the jet's propagation. Matter continue to stream into 
the cocoon only from the forward shock, and it generates a relativistic 
pressure in the cocoon. 

• The value of L can change with time, even for a constant Lj and 6*0, ac- 
cording to the slope of the density profile. In a density profile pa ~ z~"' 
with a > 2, L increases with time. This corresponds to an acceleration 
of the jet's head. In such cases a collimated jet opens up and become 
un-collimated once L becomes larger than ^q" . The opposite evolu- 
tion take place when a < 2. In the case of a = 2, L is constant, the 
head maintains its velocity, and the opening angle of a collimated jet 
remains constant. 

• The choice of the initial injection radius of the jet, z^, can affect the jet's 
breakout time, where larger values of z^ increase this time. This effect 
becomes important when the ratio of 2;* to the stellar radius > 10^'^. 
Numerical simulations typically use values of z^ which are above this 
limit and therefore their resultant breakout times are affected by this 
choice. 

Our model provides a general frame to study the properties of relativistic 
jets in different media. It can be used to examine various phenomena, such as 
the minimal energy required by the jet to break out of a boundary surface at 
a finite distance, like the edge of a star. It can also be used to test the energy 
feedback into the stellar envelope in the case of a GRB jet that penetrates 
a star, or the IGM in the case of AGN jets. This may help understands 
better issues such as the problem of IGM heating. Our model confirms that 
GRB jets, which form inside a star, are collimated before they break out 
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and become uncoUimated afterwards. Microquasars and X-ray binary jets, 
however, may start uncollimated, but the interaction with the ISM leads 
to their coUimation beyond ~ 2 ■ 10""^ pc. We also show that quasar jets 
are coUimated as well at large distances from their sources (> 10 kpc). We 
stress that in our calculations we assume that the magnetic fields in the 
jet, in the region where the jet undergoes the collimation and above, are 
dynamically unimportant and therefore can be ignored. High magnetization 
in this region will alter our results, for example by preventing the formation 
of the collimation shock, but will keep the other basic properties of the model, 
i.e. the formation of the cocoon and the collimation of the jet. 
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simulations. We also thank Y. Lubarsky for useful discussions, and to the 
anonymous referee for helpful comments. O. B. and T. P. were supported by 
the Israel Center for Excellence for High Energy Astrophysics and by an ERG 
advanced research grant. E. N. was supported in part by the Israel Science 
Foundation (grant No. 174/08) and by an EU International Reintegration 
Grant. R. S. is partially supported by IRG & ERG grant, and a Packard 
fellowship. 

APPENDIX A - Testing the assumption of en- 
ergy conservation in the jet 

In calculating the parameter L and the velocity of the head, we assume that 
the jet material does not lose energy to work as it flows from the injection 
point to the head. Such an assumption is natural if the jet is conical, since 
the jet has a constant opening angle and its envelope at a given point doesn't 
expand with time. However when the jet is coUimated its width changes with 
time, and the jet loses or gains energy due to mechanical work preformed 
against the pressure of the cocoon. When L -C 1 the head moves at a sub- 
relativistic velocity. In this limit, to a first order in /3/i, the work done due to 
sideways expansion of the jet is dW = PcdV = PcZhd'Ej. This work can be 
neglected as long as it is much smaller than the energy added to the system: 
Ljdt. We therefore define e^ = -f^-^ as the relative amount of energy 
lost to expansion of the jet. Our approximation holds as long as ej <^ 1. 
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Substituting Sj from eq. (TT2|) and Pc from eq. ( TT6|) gives: 

6, ^ ^o'/?/. = 0.03 (^^) /?,. (A-1) 

Implying that the assumption of neghgible energy losses in the jet is always 
valid in the limit of small injection angle. 

When 1 <^ L <^ 9q the head is relativistic, and the energy flux that 
enters the head through the reverse shock is reduced by a factor of ~ F^^. 
Since we are interested in the energy that enters the head we compare the 
work lost to expansion with LjTj^^dt. This gives e^- = ^^^^ -^i and applying 
the same consideration as before we get that 

6, ^ Tlel (A-2) 

In this regime Th < 0^ (see table 1), therefore the relative energy lost to 

4/3 

work in a relativistic coUimated jet maintains e^ < 6q . But in this regime 
of a collimated jet with a relativistic head, the injection angle is limited by 
Oq -C 1. This guarantees that e^ <ti 1, ar 
is always valid in this type of jets as well. 



Oq -C 1. This guarantees that e^ ^ 1, and implies that our approximation 



APPENDIX B - Analytic solutions to the rel- 
ativistic and non-relativistic limits 

The system's behavior is determined by the equations (121151161 fT2l) . Generally 
these equations should be solved numerically, due to the non-linearity of /?/,. 
But in the limits of L ^ 1 and L ^ 1, eq. ([2]) can be linearized and the 
integration over time can be solved analytically. We define the following 
integration parameters: pa{zh) = f PadV/V = gpa{zh), z^ = J f3hdt = C/3ht, 
Tc = J Pcdt = iPct and and / ^h^dt = eTj^'^t. Using these parameters, eq. 
([7]) can be written as: 



Each of these parameters takes a different value when L ^ 1 and when 
L ^ 1. If the density profile is a powerlaw of the sort pa ~ 2;~", the 
parameters become constants and their value depends on a. The resulting 
solutions including the integration parameters in each regime are presented 
below, for convenience in presenting g we took density profiles with a < 3. 
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A collimated jet with a non-relativistic head (L <C 1) 



In this regime ^ = 1 and the rest of the integration parameters takes the 



following values: C = ^ 
parameters is: 

Zh -- 



^'^ 



(3-q) • 
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A coUimated jet with a relativistic head (1 <C L < ^q ) 

In this regime ( = 1, ^ = y^^, e = g^, and g = tt^- The solutions to the 
system's parameters is: 
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Zh ~ ct 
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(B-12) 
(B-13) 

(B-14) 

(B-15) 

(B-16) 

(B-17) 
(B-18) 

(B-19) 

(B-20) 
(B-21) 
(B-22) 



The jet is collimated as long as 5 < Zh. Substituting eqn. (B-12, B-15) 
in eq. (B-22), and using the relation L = ^r^ (eq. [E]), we get that the jet 
is collimated as long as 



L<a 



-4/3 







167r3/2^e 



3e^ 



2/3 



(B-23) 



An un-collimated jet with a causally connected relativis- 
tic head {6^^^^ < L <C 6^^) 

The integration parameters in this regime are the same as in the relativistic, 
collimated regime, i.e. ( = 1, ^ = t^, e = ^^ and , g = t^A-t- But since 



(3-«) • 
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the jet is now conical the solutions are: 



Zh ~ ct (B-24) 

Ph ^ 1 (B-25) 

(3c = ( ^y ( At-Y' c-'^' (B-28) 



^' ^ l^j (lei) ^' '^-'^' 

^c = Pc (B-30) 

^j- = ^0- (B-31) 

An un-coUimated jet with an un-causal relativistic head 

(4^0^' « L) 

In this regime different parts in the head are not in causal connection and 
therefore only a fraction 2/(r/j6'o) of the energy in the head goes into the 
cocoon. The cocoon's pressure is relativistic and it pushes the edge of the 
cocoon to a velocity /3c — !■ 1 which implies that our approximation of the 
cylindrical cocoon breaks down. Moreover since the (3c > c/vS, the cocoon 
is no longer causally connected in the lateral direction and the pressure is 
no longer uniform. In this case we are able to provide the average energy 
density in the cocoon, by taking the total energy that enters the cocoon and 
dividing it with Vc, the cocoon's volume approximated as a sphere in this 
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regime. Under these approximations the model parameters are: 

(B-32) 
(B-33) 

(B-34) 

(B-35) 

(B-36) 
(B-37) 
(B-38) 

When L ^ 4rjo the Lorentz factor of the head is ~ Fjo, and the reverse 
shock becomes Newtonian. In this hmit most of the energy that flows into 
the cocoon comes from the material behind the forward shock which remains 
relativistic. Energy enters the head through the forward shock at a rate: 
Eah = ^^^^or^oPaC^5 ^^d ^ fraction 2/{Vjq6q) of that flows into the cocoon. 
Approximating the cocoon's volume as a sphere we get: 

Th = r,o (B-39) 

Y = l^JoOoPac', (B-40) 

where all the rest of the parameters remains as before. Note that in this 
regime the average energy density in the cocoon remains constant. 
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